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ABSTRACT 


An  adaptive  control  scheme  may  provide  the  best  approach  to 
the  problem  of  accommodating  automatic  flight  control  systems  to 
the  variations  of  dynamic  characteristics  encountered  over  the  flight 
envelope  of  the  aircraft.      The   C*-  Criterion  for  aircraft  time   re- 
sponse provides  a  basis  for  the  design  of  such  a  system. 

To  provide  a  margin  of  safety,    all  control  systems  have 
redundant  channels  for  emergency  use.     An  adaptive  control  system 
can  be  designed  to  be  self-organizing,    and  in  such  a  configuration 
can  provide  its  own  failure  monitoring,    resulting  in  a  simpler  and 
more  efficient  system, 

This  study  shows  that  such  a  system  is  feasible,  that  the 
response  of  the  system  is  within  the  limits  set  forth  in  the  C*- 
Criterion,  and  that  the  self-organizing  characteristics  provide 
reliable  operation  over  the  whole   range  of  flight  operations. 
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I.       INTRODUCTION 

A  self-organizing  system  may  be  defined  as  one  which  changes 
its  basic  structure  as  a  function  of  its  experiences  and/or  environ- 
ment,   with  the  general  aim  being  to  evolve  toward  some  desired 
state  of  output  behavior  or  mode  of  operation.     Such  a  system  can 
be  described  as  having  the  ability  to  learn,    and  can  therefore  be 
designed  to  handle   situations  which  rigidly  programmed  systems 
are  not  particularly  adept  in  handling.      A  report  by  Gibson,    Fu, 
et  al.  ,    Ref.    1,    gives  a  good  introduction  to,    and  provides  back- 
ground information  on,    learning  control  systems.     A  further  and 
more  recent  introduction  to  the  subject  may  be  obtained  from  a 
report  by  Mendel,    Ref.    2. 

Almost  all  of  our  present  aircraft  use  some  type  of  air-data 
scheduling  as  the  method  of  accommodating  their  automatic  flight 
control  systems  to  variations  in  the  dynamic  characteristics  en- 
countered over  the  full  flight  envelope  of  the  aircraft,      Such 
systems   rely  on  measurements  of  dynamic  pressure  or  Mach 
number.     Also,    each  system  must  be  tailored  for  the  particular 
aircraft  in  which  it  will  be  deployed,    with  resulting  extensive  design 
changes  from  system  to  system,    and  expensive  flight  testing  for 
final  adjustments  to  insure  adequate  performance. 

With  the  rapid  increase  in  technology  in  the  field  of  electronics, 
many  new  electronic  instruments  have  appeared  and  proven  to  be 
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almost  indispensable  in  the   safe  operation  of  aircraft.      These  in- 
struments,   already  available  and  in  use,    can  provide  the  necessary 
measurements  for  an  adaptive  flight  control  system.      Such  a  system 
would  depend  only  on  measurements  of  dynamic  performance  to 
provide  the  necessary  gain-programming.      This  should  give  im- 
proved system  performance  as  compared  to  a  programmed- gain 
system.      Since  the  self-adaptive  technique  uses  a  model  reference, 
it  should  be  possible  to  make  such  a  system  almost  universally 
applicable,    requiring  only  a  change  in  the  model  and  a  few  parameter 
adjustments  for  each  different  aircraft. 

The  method  of  approach  used  in  this   study  is  to  extend  that 
used  for  the  North  American  SIDAC  controller  developed  by 
Shipley  and  his  associates,    Refs   3  and  4.      This  was  modified  by 
Rang,    Ref.    5,    from  the  observation  that  the  equation  for  the 
handling  qualities  given  by  the   G*- Criterion,    which  was  developed 
by  Tobie,    Elliot,    and  Malcom,    Ref.    6,    is  very  similar  to  the  basic 
short  period  equation  of  motion  of  the  aircraft.     A  feedback  system 
using  variable  gains  and  fixed  values  representative  of  the  ranges 
of  the  aircraft  coefficients  over  the  flight  envelope  were  used  to 
meet  the   G*-  Criterion  requirements,      The  gain  changing  mechanism 
is  found  by  a  gradient  technique.      The  system  is  then  recast  into 
one  which  should  show  self-organizing  properties.      The  development 
of  the  following  theory  is  by  Shipley  [Ref.    4].      It  is  of  the  basic  form 

and  shows  the  essential  ideas.      The  extension  to  a  self-organizing 
configuration  is  given  in  the  Appendix. 
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II.       THEORY 


Beginning  with  the  short  period  perturbation  equations  for 
straight  and  level  flight,    which  may  be  written  as 

(i)     e    =    /n,  e    +    m^  «    +  ^  a    +   m^  s«. 
(2)      n    *    l/0  (  e  -  <*)     x      -  z^  <x    -  zs  st 

angle  of  attack  is  eliminated  because  it  is  difficult  to  measure 
accurately  and  its  measurement  is  affected  by  gusts  and  air  tur- 


bulence. 
(3)      0 


Equations  (1)  and  (2)  may  then  be  written  as 

j-%;  -  H^;<\    +    Mfc  [e-  %. 


m,  e    +     m* 


(4)     n   --     -^  fe  -  %.]      -    2LS  i. 


st  -e 


The  coefficients  are  abbreviated  as 

These  coefficients  are  negative  for  all  flight  conditions.      Their 
values  are  listed  in  Table  I.      Equation  (3)  then  becomes 


(5)     -  ^  e    +  41  e 


n.+  i** 


=      0 


The  C*-  Criterion  for  aircraft  response  is  given  by 


(6)  C 


*      _ 


y\ 


XP  e      *    \j<  e 


This  criterion  defines  the  gravity  forces  that  a  pilot  is  subjected 
to  as  the  aircraft  responds  to  control  inputs. 
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If  C*  is  required  to  be  a  multiple  of  the  command  input  C    ,    then 

(7)        C*       "-       KCC 

and  the  criterion  can  be  written  as 

(8)       ip  e     +    uc  e     -     n     -    K  ct      -    o 

Equations  (5)  and  (8)  are  of  the  same  form  and  are  exact  if  their 
coefficients  are  made  to  be  equal.      Thus 


To  make  the  problem  as   simple  as  possible,    the  effects  of 
actuators,    of  any  feedback  through  the  aircraft  control  system,    and 
of  the  dynamics  of  the  control  system  were  neglected.      This  ideal 
system  is  diagrammed  in  Figure   1.      It  was  then  assumed  that  the 
elevator  command  was  a  linear  addition  of  the  stick  command  and 
the  feedback.      Then 

(9)  oe       =        C       *      f 

where  the  feedback  was  a  function  of  the  gains,    so  that 

(10)  { :    «       V,  G      -     Vn  n      +     ^     C 

and  the  gains       \c     ,       \*    f    and      \g     were   such  that  there  was  no 
system  error.      The  system  error  was  defined  as 

(11)  £     =       9     -    /?,  e     -    &A  n    -     £,   C 

When  (10)  was  substituted  into  (9),    and  the  resultant  equation  and 
(5)  were  together  substituted  into  (11),    the  error  became 
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(12)        Q 


A,  -  X   -  A  Z 


s     ^   J 


e 


~#A   -A,  *  A^]   n 


In  order  to  hold  the  error  to  zero,    the  coefficients  of     ©  ,     Yl     , 
and  C  must  be  zero.      Therefore 


V. 


a. 


Z  -- 


/3n   -  /3„ 


/** 


,   vs- 


Values  of  these  coefficients  with 


4   -      -3 


^n  "-     -^       i       &x 


-%Q 


are  listed  in  Table  II. 

The  feedback  system  was  composed  of  three  separate  feed- 
back calculators,    each  set  for  a  given  flight  condition.      These 
feedbacks  were  multiplied  by  gains  K,    and  combined  to  give  the 
total  system  feedback.      The  method  of  calculating  these  gains  was 
developed  by  Shipley,    [Ref.    4]. 
Letting 

(13,  XV  =  ^  B,1    *   t  B„*     ♦    i   B; 


where 
8 


1 

Bn      = 


A,      ~    *, 


^s    V, 


^    -  A     ^  ^s  V„ 


Bs    =     A    "  A 


^s   ^ 


and      Re,     ,      K^    ,    and      Ks     are  positive  numbers, 
Then  differentiation  with  respect  to  time  gave 
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<14>  ft 


-~      A 


A  B,#    +i  Bn£    -  £  0 


g     J1- 


and  choosing 


7T1     =       "i   ©  & 
d^w      =        k    vn    G 

Jt 


K<   C    G 


resulted  in  the  equation 
(15)  ft  =  /?s  6  G 
Since    ]5%  <  0  t       G  was  chosen  so  that 


C-     -        <.«f\    € 


4  /  ,     0  £  e0  <  e 

0  ,    -€0  ^  e  <  £0 

-I  £    <  -  £0  £0 


which  insured  that        -pp         was  negative.     Therefore    V    will 
always  decrease  toward  zero  with  time  and  the  gains  tended  toward 
an  ideal  value  for  each  flight  condition.      In  this   system  the  gains 
are  calculated  from  the  equations 


with 


-5 


S      .     I<„    =     5"  X  10".        K,    =    S'OO    .       €£ 


5  M 


oz 
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III.       EXPERIMENTAL  PROCEDURES 

The   system  was   set  up  as  shown  in  Figures  2  and  3,    using 
the   COMCOR  5000  analog  computer.      The  basic  aircraft  analog 
was  tested  by  making  plots  of   0    and   vn    as  a  step  input    S^    was 
applied  to  the  system.      The  values  obtained  were  compared  with 
computed  values  of  steady  state  amplitude,    frequency,    and  damping 
ratio.      Transfer  functions  were  derived  for    @   and    ^     with  respect 
to  a  step  input,    and  the  final  value  theorem  was  applied  to  obtain 
the  steady  state  ratios. 

After  the  basic  aircraft  analog  was  tested,    the  feedback  system 
was  incorporated.     Each  of  the  three  feedback  elements  was  chosen 
to  represent  a  particular  flight  condition.      The  three  flight  con- 
ditions chosen  were  0002,    0009,    and  5020,    where  the  first  two 
numbers  in  each  series   represent  the  altitude  in  thousands  of  feet 
and  the  last  two  represent  the   Mach  number  in  tenths.      By  choosing 
these  particular  flight  conditions  the  operating  range  of  the  air- 
craft was  fairly  well  represented.      Flight  condition  0002  is  the 
low  altitude,    low  airspeed,    landing  condition.      Flight  condition 
0009  is  a  subsonic,    low  altitude  condition  where  the  aircraft  has 
excellent  response  to  control  inputs  and  a  high  damping  ratio. 
Flight  condition  5020  is  a  high  altitude,    high  speed  condition  where 
the  control  response  is  not  very  good  and  the  damping  ratio  is 
low. 
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The  three  feedback  elements  were  tested  by  setting  the  value 
of  K  for  the  element  being  tested  at  1.  0,    and  the  K's  for  the  other 
two  elements  at  zero.      Since  the  element  was  set  up  to  represent 
exactly  the  flight  condition  being  considered,    the  values  of  K 
should  not  change  and  the  error  should  be  zero.      This  proved  to 
be  the  case  for  each  of  the  three  basic  flight  conditions  when  K 
was  set  in  as  an  initial  value.      However,    when  operating  the  system 
with  no  initial  values  set  in,    the  K's   reached  a  steady  state  value 
which  was  somewhat  different  from  the  combination  of  a  1.0  and 
two  zeros.      Each  flight  condition  did  reach  a  steady  state  in  which 
the  error  signal  remained  in  its  dead  band.      For  an  example,    see 
Figures  10,    11,    and  12. 

After  both  the  aircraft  analog  and  the  feedback  analog  were 
tested,    the  system  was  tested  for  self-organizing  characteristics. 
For  simulation  purposes,    the  values  of    0  ,     W    ,    and  C  were  all 
taken  from  the  aircraft  analog.      In  a  real  system,    these  values 
would  come  from  separate  sensors.     Sensor  failures  were  simu- 
lated by  changing  the  value  of      \  =     ,      \n      ,    or     \s     in  one  of  the 
feedback  elements  from  its  calculated  value  to  a  different  value. 
In  this   simulation,    the  values  for  errors  were  zero  and  one. 

(id  f    --  X^  e   •*  V.  r\  -  Vs   c 

Setting  the  value  to  zero  simulated  the  loss  of  signal  input  from 
a  sensor,  while  setting  the  value  to  one  simulated  an  erroneous 
signal  input.      One  other  type  of  error  which  was  simulated  consisted 
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of  using  a  steady  input  in  place  of  one  of  the  variable  signals  for 
Q     ,     n     ,    and  C  from  the  aircraft  analog.     This  would  simulate 
a  sensor  which  was  stuck  and  producing  a  constant  output,    regard- 
less of  the  actual  conditions.     This  last  type  error  was  not  tested 
extensively,    since  the  proper  values  for  such  signals  were  not 
known. 
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IV.       RESULTS 


The  full  system,    for  all  flight  conditions  tested,    drove  its 
gains  to  values   such  that  the  system  error  remained  in  its  zero 
range.      The  time  required  for  this  system  adjustment  varied  with 
flight  condition,    frequency  and  amplitude  of  the  input  commands, 
and  the  number  of  commands  given.      The  flight  characteristics 
did  not  vary    significantly  during  this  transient  time.     (See  Figures 
13,    14,    and  15) 

The  self-organizing  characteristics  of  the   system  appear 
to  be  adequate.     Simulated  failures,    as  outlined  in  the  method  of 
testing,    were  compensated  by  variations  in  the  gains.      The  plots 
in  Figures   16  thru  21   show  little  change  in  aircraft  flight  character- 
istics.     In  some  cases,    the  value  of  the  error  signal  could  not  be 
made  to  remain  in  its  dead  band,    but  always  quickly  drove  toward 
zero.      For  the  three  design  flight  conditions,    any  simulated 
failures  outside  of  the  design  feedback  element  showed  no  effect 
on  flight  characteristics.      This  is  as  expected,    since  the  design 
element  should  handle  the  flight  condition  by  itself.      When  the 
failure  was   simulated  in  the  design  feedback  element,    there  was 
usually  a  big  change  in  the  variable  gains,    and  a  slight  change  in 
performance.      The  flight  characteristics  were  still  satisfactory. 

When  multiple  failures  of  the  zero  input  type  were  simulated, 
the  system  was   still  able  to  function  adequately.      There  was  some 
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deterioration  in  the  amplitude  of  response  and  the  time  response 
characteristics.      The  system  error  usually  did  not  stay  in  its  zero 
range  'when  command  signals  were  given,    but  always  went  rapidly 
to  zero.      The  system  was  able  to  handle  combinations  of  up  to  six 
simulated  failures  of  inputs,    as  long  as  at  least  one  each  of  the    0  , 
Y\      ,    and  C  inputs  remained. 

When  multiple  failures  were  simulated  using  erroneous  input 
signals,    system  reaction  was  not  as  good  as  with  the  zero  input 
failure  simulation.     The  system  was  still  able  to  reach  a  semi- 
steady  state  condition  in  which  the  error  would  drive  to  zero,    but 
response  was  generally  poorer.     The  poorest  results  were  obtained 
when  the  error  simulation  involved  two  of  the  same  quantities, 
such  as  two     Q  s.     In  such  cases,    the  system  could  not  tell  which 
input  was  the  correct  one,    and  the  resulting  performance  plots 
were  quite  bad.     Similar  results  were  obtained  by  mixing  the 
modes  of  failure  together.      It  is  considered  that  part  of  the  problem 
encountered  was  the  result  of  testing  techniques,    in  which  the 
inputs  were  increased  by  as  much  as  a  factor  of  seven  to  provide 
the  erroneous  signals. 

One  other  method  was  used  to  simulate  failures.     A  constant 
signal  was  used  to  replace  one  of  the  inputs.     The  system  was  able 
to  handle  this  type  of  error  very  well.      This  type  of  failure  was 
not  tested  extensively  since  representative  error  values  were  not 
known. 
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V.        CONCLUSIONS 

A  self-organizing  adaptive  control  scheme  can  handle  gain 
scheduling  within  a  control  system  in  such  a  way  that  the   system 
response  is  relatively  invariant  as  required  by  the   C*-  Criterion. 
The  ability  to  alter  system  gains  as  necessary  to  meet  the  response 
criterion  gives  the  system  better  response  than  can  be  realized 
with  a  fixed  gain  schedule  system,      The  system  is  also  able  to 
accommodate  sensor  failures,    thus  providing  a  built-in  backup. 
Since  the  system  uses  a  model  reference,    it  should  be  easily 
adaptable  to  all  types  of  aircraft  by  changing  the  model  and  the 
fixed  feedback  parameters  to  those  of  the  particular  aircraft. 
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APPENDIX 


The  equations  of  motion  of  the   system  are 


(i) 


I 


K 


e  =    4,  e   +  ^  n   +    ^8  se 

Sfc  =  C  +      -f 


Changing  the  configuration  from  one  with  a  single  feedback     f    , 


which  has  three  varying  gains      *\"    ,   VM    ,      Vs 
has  three  feedbacks  of  the  form  of      T    , 


,    to  one  which 


(2) 


but  with  fixed  gains.      In  matrix  form 


(3) 


V 


r; 

v; 

c 

C 

V3 

v: 

v; 

is  constant.     The  three  feedbacks  are  combined  with  varying 


gains. 


(4) 


K    = 


K. 


K 


K 


to  give  the  final  form  of  the  input  into  the  elevator  channel 
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(5)        -f     =       K.-f,      +    Kxfx        +    K,f, 


To  facilitate  the  derivation,    change  to  vector-matrix  notation. 


Let 


X     = 


e 

Y) 
C 


/3     = 


S   J 


z?    - 


/3S 


*   = 


The  set  of  equations  from  (1)  then  yields 


(6) 


where  the  prime  denotes  the  transpose  of  the  matrix. 


Consider  the  function 


(7)      2\/   -(*'-*'+  3*  *")Q  (■*-*  +  A  *"  )    . 


where  Q  is  a  positive  definite  matrix  to  be  chosen  later, 
Differentiation  gives 


(8,        £     - 


7*    jy 


(£'-*'  +  &  f)  a  zs  # 


Since     V    is  to  be  driven  to  zero,    which  in  turn  means      €      is 


driven  to  zero,    the  requirement  is 


IFF 


£    0 


Thus,    choose 


(9) 

so  that 


4J  = 


_i 


Q     x  G 


/S,    £    G- 
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Since  the  parameter      P%      is  always  negative,    taking 


G     =         *$ 


5Gn      € 


makes 


dt 


£    |i| 


so  that 


If 


&  «  ° 


(10) 


Q' 


K,         O        O 
O  K*       O 

O         O       J<$ 


then  the  system  is  the  same  as  given  in  the  body  of  the  thesis. 
Equation  (5)  combined  with  (2)  is 


(11)  f     =        !<'  V   x 


But  also  from  equation  (1) 


(12) 


y    x 


Therefore 


r  =     V'k 


Equation  (9)  is  then 


QT'x  G 


or 


(13) 


(QTT'  x   G 
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Now  choose 


(  Q  V) 


-i 


Qll  V 


so  that 


(14) 


with 


JK 
Jt 


q: 


a 


-i 


^     J 


G 


K1 

o 

o 

o 

k. 

0 

0 

o 

K*_ 

Equation  (14)  becomes 


f  ^-     -_         v     f,  G 
Jt-  ^ 


(15) 


.  Jt 


These  produce  the  variable  gains  used  in  the  self-organizing 


system.     Since 


Q 


V  Qo  fv') 


i  ^\  ' 


If      \      is  non-singular  and    Q0       is  positive  definite,    then  Q  is 
positive  definite,    which  is  necessary  for  concluding  that      €.       goes 
to  zero  as      V    goes  to  zero. 
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FLIGHT 
CONDITION 

0002 

0009 

1504 

3006 

4509 

5020 

ALTITUDE 

0 

0 

15,000 

30,000 

45,000 

50,000 

MACH  NO. 

0.2 

0.9 

0.4 

0.6 

0.9 

2.0 

Uo 

223 

1005 

423 

597 

872 

1938 

-nm 

1.3 

38 

2.93 

3.6 

4.35 

42.4 

-   Mi 

.26 

1.38 

.33 

.27 

.20 

.04 

-  M, 

.43 

2.61 

.56 

.51 

.34 

.48 

-  ** 

2.8 

58.4 

6.3 

7.5 

6.52 

14.9 

-  2« 

83 

2420 

228 

262 

277 

817 

-  2. 

13.2 

313 

31.5 

38.2 

35.8 

102 

-  ^ 

2.58 

53.5 

5.90 

6.93 

3.78 

9.61 

"  *, 

.69 

4.0 

.89 

.78 

.54 

.52 

-    A, 

.0146 

.0143 

.0121 

.0133 

.0158 

.0519 

^, 

90 

-.019 

.36 

.32 

.65 

.26 

^ 

.0021 

.0001 

.0013 

.0010 

.0011 

.0033 

^ 

6.75 

-.63 

2.39 

1.88 

4.29 

1.08 

TABLE  I  AIRCRAFT  PARAMETERS 

(F4  AIRCRAFT  Units  in  radians,  feet,  seconds) 
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K, 

,05 

Kn 

.0005 

"<s 

500 

,895 

< 

-.0187 

^5 

.258 

v: 

.00213 

v- 

.000107 

v: 

-.00331 

v; 

6,75 

c 

-.626 

vs3 

1.083 

£c 

.02 

*1 

-3.0 

A, 

-.02 

£6 

-20 

TABLE  II 


FEEDBACK  CONSTANTS 
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CM 
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FIGURE  4 

Plots  for  flight  condition  0002,    showing  the  response  of 
the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  5 

Plots  for  flight  condition  0009,    showing  the  response  of 
the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  6 

Plots  for  flight  condition  5020,    showing  the  response 
of  the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  7 

Plots  for  flight  condition  1504,    showing  the  response 
of  the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  8 

Plots  for  flight  condition  3006,    showing  the   response 
of  the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  9 

Plots  for  flight  condition  4509,    showing  the  response 
of  the  basic  aircraft  alone  and  with  feedback. 
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FIGURE  10 

Plot  for  flight  condition  0002,    showing  K's  driving 
toward  their  steady  state  values. 
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FIGURE  11 

Plot  for  flight  condition  0009,    showing  K's  driving 
toward  their  steady  state  values. 
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FIGURE  12 

Plot  for  flight  condition  5020,    showing  K's  driving 
toward  their  steady  state  values. 
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FIGURE  13 

Plot  for  flight  condition  0002,  .showing  the  change  in  response 
as  the  K's  go  from  zero  toward  their  steady  state  value. 
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FIGURE  14 

Plot  for  flight  condition  0009,    showing  the  change  in  response 
as  the  K's  go  from  zero  toward  their  steady  state  value. 
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FIGURE  15 

Plot  for  flight  condition  5020,    showing  the  change  in  response 
as  the  K's  go  from  zero  toward  their  steady  state  value. 
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FIGURE  16 


Plots  for  flight  condition  0002,    with  the  value  of      VV     .     V„ 

r->3  7    ,  ~ 

and      \1    .respectively,    set  to  zero. 
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FIGURE  17 


Plots  for  flight  condition  0002,   with  the  value  of     V„     ,     V^     , 
and      \  n     ,  respectively,    set  to  zero. 
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FIGURE  18 


Plots  for  flight  condition  0002,   with  the  value  of     YV  V? 

i     ,    respectively,    set  to  zero. 
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FIGURE  19 


Plots  for  flight  condition  0002,   with  the  value  of     Y),      t     V,      , 
and      \       »    respectively,    set  to  one. 
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FIGURE  20 


Plots  for  flight  condition  0002,   with  the  value  of    V„     ,     Y„ 
and      V„      ,    respectively,    set  to  one. 


48 


-/.<]- 


9 


■»!.« 


-1- 


♦  .1 


■1.0 


W.d 


■10 


/» 


-'£) 


i        i        i      ~r 1 1 i       t 1    | p 


2: 


5 


=&: 


b- 


t2 


-+  I — l — I- 


f.1 


t 


H 1 1- 


-I 1 (- 


f^ftfcr 


\p- 


^ 


£ 


-| 1  I 1 1 1 1 1 1 1 1 1 


H 1 h 


H *  I 1 1 1 1 1 (- 


$ 


£ 


j — i — i — i — i — i i u_ 


I 1 1 1 1 1 u 


b-- 


V- 


£^ 


77 


H 1 1 \ — I 1- 


*F 


H 1 1 h 


^ 1 h 


FIGURE  21 


Plots  for  flight  condition  00.02,    with  the  value  of    V4    ,      Vj1     , 
and      \  r     ,    respectively,    set  to  one. 
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